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Abstract 

We calculate the effective potential of A^ 4 theory at finite tempera- 
ture under the super-daisy approximation, after expressing its deriva- 
tive with respect to mass square in terms of the full propagator. This 
expression becomes the self-consistent equation for the derivative of 
the effective potential. We find the phase transition is first order with 
this approximation. We compare our result with others. 
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1 introduction 



The symmetry restoration, or the phase transition, at high temperature is re- 
garded today to play an important role in particle physics and cosmology [Q]. 
The investigation of the phase transition, however, is found difficult quite 
often because of the unreliability of the perturbation theory at high temper- 
ature. In case of the electroweak phase transition, which may be important 
for the baryogenesis@, for example, such a breakdown of the ordinary per- 
turbative expansion occurs for the Higgs boson mass m# to be equal to or 
larger than the weak boson mass M^[0. The experimental evidences 0, in 
fact, suggest that uih is not smaller than Mw, indicating that we cannot 
examine the electroweak phase transition by the perturbation theory. 

The difficulty is caused by the Bose condensation in the thermal bath, and 
is a well-known feature inherent to the finite temperature field theory^, 
The traditional method to improve the perturbation theory is to resum the 
daisy diagrams^ [7|, |L [J |TI], [TT], [T2|, |Tj| and use the mass m?, the sum of the 
thermal mass and the zero temperature mass, instead of the latter only. But 
this improvement is not sufficient. When one uses the loop expansion, one can 
estimate how small a contribution the next order loop has in comparison with 
the previous order loop|| ||. The ratio of these is called a loop expansion 
parameter. After resumming the daisy diagrams, we find the loop expansion 
parameter is in the A(/> 4 theory for the temperature T. It is 0(1) near the 
critical temperature and the analysis with the daisy diagram indicates the 
phase transition to be of first order though it is well known to be of second 
order 0, [TJ|]. This shows the perturbation theory breaks down near the critical 



temperature and is not reliable. Many non-purturbative approaches, such 
as lattice Monte Carlo [Io| , epsilon expansion flIo|, CJT method WJ^ |T8| , |T5| , 



effective three dimensional theory pOfl, gap equation met ho d [pl}1 etc., have 



been made to investigate the phase transition at high temperature. However, 
they have not reached a consensus on what parameter region of and 
Motile phase transition is of first order. 

Recently, two of the present authors(K.O. and J.S.) proposed the new 
non-perturbative method to calculate the effective potential V^[p2] . We ex- 



press its derivative with respect to the mass square, dV/dm 2 , in term of 
the full propagator. We calculate the effective potential by integrating the 
derivative of the effective potential with an initial condition at m 2 = m 2 , 
given in the region where the perturbation theory is reliable: 

V(m 2 = -/i 2 ) = V{m 2 = m 2 ) + / Tr^dm 2 , (1) 



dm 2 

where — fi 2 is the mass parameter of the theory and m is the mass scale 
at which we give the initial condition V(m 2 = r/ig). In principle we can 
get the effective potential non-perturbatively by this procedure. The main 
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problem of this method is how to approximate the full propagator in it. In 
the previous paper we ignored the momentum dependence of the self energy 
graphs and replaced it by the second derivative of the effective potential with 
respect to the expectation value of the field. In this paper we try another 
approximation of the full propagator. We can sum up all the super-daisy 
diagram contributions correctly in this methodP]. We evaluate the effective 
potential numerically with this approximation. 

The paper is organised as follows. In Sec. 2 we discuss the structure of 
the evolution equation that we have to solve. In Sec. 3 we introduce the 
new approximation. In Sec. 4 we show the numerical result. In Sec. 5 we 
summarise and discuss our result, comparing our method with the others. 



2 The structure of the evolution equation 

In this section we discuss the evolution of the effective potential with respect 
to the mass. We consider the A0 4 theory which is defined by the Lagrangian 
density 



1 

2 



dr 



l -{V<Pf+ l -^-^ + J<P + c.t. (2) 



In the previous paper |22[ we derived the evolution equation of the effective 
potential, 
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where II = II(p 2 , — Pq, <p, m 2 , r) is the sum of all the one particle irreducible (1PI) 
self energy graphs. 

The first term of Eq. (|3]) is the simple background field contribution. 
The second term gives the finite temperature contribution. The third term 
expresses the effect which remains finite at zero temperature. The last term 
is the counter term contribution. 



1 One can also do this using CJT method |]l7|. We will mention this in SecJ|. 
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Figure 1: 1 PI self energy graphs of II. The line with striped blob repre- 
sents the full propagator. The circle with a net expresses the full 4-point 
vertex(the third term in RHS). The circle with bricks represents the full 3- 
point vertex(the forth term in RHS). The last with a black blob term of RHS 
represents the counter term. 
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(a) (b) 

Figure 3: Examples of super-daisy diagrams. The graph (a) is often called 
merely a daisy diagram. 

3 Super-Daisy approximation 

The main problem of our method is how to approximate II in the Eq. (|3]) . The 
function II is expressed by the full propagator, the full three and the full four 
point functions (Fig.|]) according to the Schwinger - Dyson equation. The 
Schwinger - Dyson equation expresses the inverse of the full propagator by the 
relation expressed by Fig.^. In the right hand side of Fig.|l|, the first and the 
forth terms are independent of the external momentum and the second and 
the third terms depend on them. In this paper we assume the simplest case 
that the momentum independent term (the first term in Fig.[I|) is dominant. 
This corresponds to the well-known super-daisy approximation |7]] which sums 
up all the graphs like Fig.[| We stress that the effective potential consists of 
all the super-daisy diagrams without over counting by this approximation. 
In this approximation we see the following relation(Fig.[|), 

n = - \n (4) 

In the following we ignore the third and forth terms in Eq.(j3|), considering 
that only the first two terms are important to discuss the phase transition^] . 
Due to this approximation, we neglect the loop contribution remaining finite 
at T = 0. Using relation (£|) and integrating Eq.(H) over p [|H], [II], |24j| , we 
obtain the approximate evolution equation, 

dV 1t 2 l 1 A 00 , P 2 1 ,~ 

dm* ~ 2 + 2Wo dP I W 1 1 / 8V\ ' {) 

2 Y 2 + m2 + X fr^ exp - Wp 2 + m 2 + - 1 
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dV = 112 + 




Figure 4: The relation between the self energy II and the derivative of the 
effective potential V with respect to the mass square. See the first two term 
ofEq.@. 



4 Numerical result 

In the previous section we have obtained the evolution equation for 
under the super-daisy approximation. In this section we explain details of 
the numerical calculation and show its result. 



4.1 Details of the numerical calculation 



We can calculate the effective potential, by solving the Eq. ((5|) and integrat- 
ing it according to the Eq. ([I]) numerically. 

In the Eq. ([!]) we set the initial scale mo as large as T so that we can 
evaluate the effective potential very well by the loop expansion. We use the 



same initial condition as was used in the previous paper p2 



V(m 2 = ml) 



dpp 2 log 



The integral in the evolution Eq. (|5|) is well defined when the effective 
mass square, m 2 + A , is real and positive. Below the critical temperature, 
however, the effective mass square can be negative or complex. We needs 
analytic continuation of the integral of the Eq.@. This can be done by 



rewriting the Eq.(j|) as in the previous paper 



mi 



dV 1 T 2 r y/g + 2M) 



9m 2 2^ 47r 2 y e 2+M -l 



where M = ^m 2 + Agj . To find <9K/<9m 2 we solve Eq®. 

As is well-known, the effective potential becomes complex at small 
below the critical temperature. It indicates instability of the state and the 
imaginary part of the effective potential is interpreted to be related with a 
decay rate of the state[[SJ. The imaginary part arises from the integral of 
the evolution equation in our method (see Eq.([?D). It is natural to suppose 
that the imaginary part of the effective potential is negative to interpret it 
with decay rate. In order that the imaginary part of the effective potential 
be negative, the imaginary part of must be positive(see Eq.([TJ)). We find 
that there are two solutions for the Eq. ([F[) and that their imaginary part 
has the opposite sign. We choose the solution of with positive imaginary 
part. 



4.2 Result 

Let us explain the numerical result. For the graphs Fig.(^) ~ (|7]) we set] 2 ] the 
four-point coupling A = 1. 

We show the real part of the effective potential near the critical tempera- 
ture in Fig.|5]. We see a small barrier between the symmetric vacuum (<f) = 0) 
and the symmetry-broken vacuum (</> ^ 0), which clearly shows a first order 
phase transition to occur. 

A first order phase transition is also indicated by the behaviour of the 
field expectation value at the stable point, since Fig.|6] shows there is a finite 
jump of it. 

We get the negative imaginary part of the effective potential below the 
critical temperature and show it in Fig.0. We can see that the magnitude 
of the imaginary part increases as the field expectation value decreases. It 
expresses that the closer the field expectation value to the origin is, the less 
stable the state below the critical temperature is. 

In the above we used the initial condition at m — itlq = T. It is required 
itlq be the same order as the temperature in order that the initial condition 
can be evaluated by the perturbation. We have calculated the effective po- 
tential with other initial conditions mo = 2T and -| and do not have found 
any appreciable change of the effective potential. It ensures the consistency 
of our calculation. 

2 We varied the value of A and found no qualitative change. 
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Figure 5: Real part of the effective potential at T = /it. We see a 
barrier to indicate a first order phase transition. 
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Figure 6: The field expectation value near the critical temperature at the 
stable point. We see a finite jump of 4>//i at t — 4.93/i to indicate a first 
order phase transition. 



8 




9 



(a) (b) 

Figure 8: Two loop diagrams that contribute to the effective potential, (a) 
is daisy-like and (b) is not daisy-like. 

5 Summary and discussion 

In this paper we have given a new method to calculate the effective potential 
of A0 4 theory under the super-daisy approximation without over counting. 
We have numerically evaluated a real part and a imaginary part of the ef- 
fective potential without using high temperature expansion. The real part 
indicates a first order phase transition though the true order should be sec- 
ond. The imaginary part indicates the instability is larger for the smaller 
expectation value of the field below the critical temperature. 

Now we compare our method with other approaches. First we compare 
it with CJT method [|l7j which can also gather up the super-daisy graphs 
without over counting. It carried out by truncating the CJT expansion at 
0(A). A first order phase transition is indicated by the calculation of this 
model with this method as well as the high temperature expansion ||T9| . 



Second we compare with daisy-improved perturbation theory. In this 
method the effective potential indicates a first order phase transition at one 



loop level |n|. There are two graphs at two loop level(see Fig.||). If one in- 
cludes the contribution of Fig.|S|(a) alone|TB[, the phase transition is still of 
first order^. A second order phase transition is indicated when the contri- 
butions of both Fig||(a) and Fig.||(b) are included^. This indicates that the 
contribution of Fig.^(b), which is not daisy-like, seems to play the important 
role to the effective potential. 

Finally we compare our result with the approximation of the previous 
paper. In the previous paper we get the effective potential that indicates 
a second order phase transition. The contributions of the second and the 

3 In rcf. only Fig.||(a) was included. 

4 In ref. ]3| they calculated the effective potential including both of Fig.||(a) and Fig.^(b) . 
They did not, however, calculate the expectation value of the field nor argue on the order 
of the phase transition with their 2-loop result. We calculated the expectation value of the 
field using their effective potential and found the phase transition second to be of second 
order. 
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third terms of Fig.|T] with zero external momentum are taken into account 
by the previous approximation. The difference of the two approximations is 
whether we include the contributions of the daisy-like diagrams only or not. 
This also indicates that non-daisy-like diagrams are important for the second 
order phase transition to be derived correctly. 

According to the above comparisons, the diagrams which are not daisy 
like seem to give important contributions. We have to take into account the 
contribution of the second and the third terms in Fig.|IJ Their momentum 
dependence, especially the dependence on the small momentum, will be im- 
portant in Eq. (^) since the infrared behaviour of the self energy plays a very 
important role for a second order phase transition^]. This method with 
seems attractive because it may bring a chance to step in the region where 
the perturbation theory breaks down. 

We finally express our thanks to A. Niegawa and T. Inagaki for valuable 
discussions and communications. 
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